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Abstract — Two classes of turbo codes over high-order finite 
fields are introduced. The codes are derived from a particular 
protograph sub-ensemble of the {d^, — 2,dc = 3) low-density 
parity-check code ensemble. A first construction is derived 
as a parallel concatenation of two non-binary, time-variant 
accumulators. The second construction is based on the serial 
concatenation of a non-binary, time-variant differentiator and of 
a non-binary, time-variant accumulator, and provides a highly- 
structured flexible encoding scheme for = 2,dc = 4) 
ensemble codes. A cycle graph representation is provided. The 
proposed codes can be decoded efficiently either as low-density 
parity-check codes (via belief propagation decoding over the 
codes bipartite graph) or as turbo codes (via the forward- 
backward algorithm applied to the component codes trellis). The 
forward-backward algorithm for symbol maximum a posteriori 
decoding of the component codes is illustrated and simplified by 
means of the fast Fourier transform. The proposed codes provide 
remarkable gains (~ 1 dB) over binary low-density parity-check 
and turbo codes in the moderate-short block regimes. 



I. Introduction 

Low-density parity-check (LDPC) codes HI constructed on 
high-order finite fields JSJ-I?) show remarkable coding gains 
with respect to (w.rt.) binary LDPC/turbo codes ||5]-|l8]. The 
gain is especially visible in the moderate-to-short block length 
(k < 1000 bits) regime IS), H), where binary iteratively- 
decodable codes tend to exhibit either high error floors or non- 
negligible coding gain losses |[9l w.rt. available benchmarks 
(e.g., 1 dB with respect to the random coding bound (RCB) 
lfTOl )Fl Ultra-sparse non-binary LDPC codes based on fields 
of order q > 6A for short block lengths allow approaching the 
average performance of random codes by 0.2 dB down to 
medium-low codeword error rates (CERs). 

In this paper, we provide a novel turbo code construction, 
which leads to non-binary turbo codes based on convolutional 
codes on a finite field ¥g of order g = 2™ > 2. The proposed 
construction bridges rate- 1/3 non-binary turbo codes and reg- 
ular {dv = 2,dc = 3) non-binary LDPC codes, where dc and 
dv are the check node (CN) and variable node (VN) degrees, 
respectively. More specifically, the turbo code construction 
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follows from a specific expansion of a {d^ = 2,dc = 3) 
protograph0 The so-obtained turbo codes have performance 
close to that of their non-binary LDPC counterparts, with 
the advantage of an efficient encoding structure. Even more 
important, the proposed construction allows building families 
of rate-compatible turbo codes with flexible block size, by 
adopting combinatorial (on-the-fly) interleaver constructions. 

Convolutional turbo(-like) codes over non-binary finite 
fields/rings have been investigated previously, e.g. in [[T4|- 
1201. With respect to the past works, the main novelties 
of our construction are listed next. Most of the previous 
contributions devote attention to fields of low order, whereas 
we focus on high order fields, e.g. of order (7 = 2™ = 256. 
Moreover, our construction is based on codes with memory 
1 (in symbols), whereas in ifTTl . |[T9l codes with larger 
memories are considered. The convolutional codes adopted 
by the proposed construction are time-variant, while in llTSl . 
lfT7l - |fT9l the feed-forward / feedback polynomials are fixed. 
To our knowledge, the only non-binary iterative schemes 
adopting time-variant convolutional codes are the irregular 
repeat accumulate (IRA) codes of ll20l . However, in 1201 a 
turbo code structure is not considered explicitly. We further 
propose an alternative serial turbo code interpretation of the 
{dy = 2,dc = 3) regular protograph ensemble which allows 
to efficiently encode rate-1/2 {dy =2,dc ^ 4) regular LDPC 
codes as a serial concatenation of a non-binary time-variant 
differentiator, and interleaver, and a non-binary time variant 
accumulator The proposed construction features a nice graph 
interpretation which provides useful insights for the interleaver 
design. High rates can be obtained by suitably puncturing 
the code symbols, whereas low rates can be easily obtained 
by using the multiplicative repetition approach of ET\ . A 
discussion on how decoding can be performed either in a 
turbo-like fashion (i.e., by iterating a BCJR decoder on the 
(j-states trellises of the component codes) or via the classical 
belief propagation (BP) algorithm for LDPC codes is provided. 
In both cases, the decoder complexity growth w.rt. the field 
order q can be limited to 0{q\og2 q) by adopting fast Fourier 
transforms. 

The contribution is structured as follows. Section HI] de- 
scribes the code structure. In Section |IIT]the iterative decoding 
algorithm is discussed. Numerical results and conclusions 
follow in Sections |IV] and [V] respectively. 



^Similarly, in 1121 . 1131 a bridge between LDPC and turbo code construc- 
tions was provided in the binary context. However, in 1121 . 1131 protograph 
ensembles have not been considered. 



II. Code Structure 

Regular {d^ = 2,dc) non-binary LDPC codes are char- 
acterized by excellent iterative decoding thresholds Js). In 
this paper, we shall focus on the {dy = 2,(ic = 3) rate- 
1/3 case. The iterative decoding threshold over the binary- 
input additive white Gaussian noise (AWGN) channel for 
the {dy = 2,dc = 3) unstructured ensembles with random 
choice of the non-zero coefficients in the parity-check matrix is 
{Eb/NoY ~ -0.25 dfiHonly 0.25 dB away from the Shannon 
limit. The turbo codes described in this paper are constructed 
on a specific protograph sub-ensemble of the {dy = 2, dc = 3) 
regular ensemble, depicted in Fig. ^a). A protograph ll22l . 
||23l is a Tanner graph Il24l with a relatively small number 
of nodes. A protograph Q = {V,C,£) consists of a set of N 
variable nodes V, a set of M check nodes C, and a set of 
edges £. Each edge e^ j G £ connects a VN Vj £ V to a CN 
Cj G C. A larger derived graph can be obtained by a copy- 
and-permute procedure; the protograph is copied K times, and 
the edges of the individual replicas are permuted among the 
K replicas. A protograph can be described by a base matrix B 
of size M X N. The element bi j of B represents the number 
of edges connecting the VN Vj to the CN Ci . The base matrix 
associated with the protograph in Fig. [Tla) is given by 



B 



1 2 
1 2 



We proceed by expanding the protograph into the derived 
graph as follows. We first replace the upper-left '1' in B with 
a K X K identity matrix, I. The lower-left '1' is replaced by 
a K X K permutation matrix H. Each of the '2' entries of B 
is replaced by a K x K cyclic matrix 



1 
1 1 








. 1 
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1 
1 1 



The (2A') X {2>K) derived graph adjacency matrix has the form 

IPO 

nop 



We then construct a parity-check matrix on Fg by replacing 
each '1' in T with a suitable element of F* = Fq\{0}. The 
2K X iK full-rank non-binary parity-check matrix is given by 



H 



i p(i) 
n p(2) 



(1) 



where 

• 1 \s a K X K matrix with non-zero entries only on 
the main diagonal {pseudo-identity matrix). More specif- 
ically, ij^i — g^^'^ for j = I, otherwise ij^i = 0. 

• n is a if X if matrix with non-0 entries only for 
^'rO)j,Vj e (l,iv). In particular, ^jXj) = gf\Vj G 
(1, K). Thus, n can be described as the product between 
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Fig. 1. Rate-1/3 protograph as parallel (a) and serial (b) concatenation. 



a permutation matrix 11 and a pseudo-identity matrix 
with diagonal coefficients gl^\g{^ , . . . , gj^Li, with 7r(j) 
being the permutation function. 

Both p(i),p(2) are double diagonal matrices with form 



/I 



(^) 




1 

(^) 



/r^ 1 



(^) 



1,2. 



A. Parallel Concatenation 

The parity-check matrix ([T]i may be interpreted as the parity- 
check matrix of the parallel concatenated convolutional code 
(PCCC). The information word u of if = k/m symbols 
in Fg, q — 2"\ is input to a rate-1, memory-1 time-variant 
recursive systematic convolutional (RSC) tail-biting encoder 
(non-binary accumulator). The first set of parity symbols p^^^ 
is obtained as 



and with 



(1) 



n(l) 



Vi e [0, K - 1] (2) 
properly initialized^ Here, 



Pi ' i9i ' ifi' '= H^i?' ^nd all operations are in Fq. The second 
set of parity symbols are obtained by first permuting the 
symbols of u according to the interleaving rule i — 7r(i). 
The permuted information word u' (with u[ = U7r(i)) is then 
fed in a second rate-1, memory-1 time-variant RSC tail-biting 
encoder The second set of parity symbols p*^^' is obtained as 



P. 



(2) 



91 K + fl 'p-Ji 



f(2)J2) 



e [0, K - 1] 



(3) 



(2) (2) 

and with p_{ = Px_i- The systematic codeword is given by 
c = [u|p(^) Ip*-^'] and has length 3if symbols. The code length 
is n — 3Km bits and the code dimension is fc = Km. The 
code rate for the proposed construction is r ^ k/n ^ 1/3. 



^Throughout the paper. Ej, will denote the energy per information bit and 
Nf) the one-sided noise power spectral density. 



^The initialization of (and consequently of can be obtained as 

indicated in 1251 Lemma 4]. 



B. Serial Concatenation 

By stretching the protograph of Fig. [Tla) onto the one 
of Fig. [lib), it is possible to devise an ahernative encoder 
which is based on the serial concatenation of a memory- 1 non- 
recursive encoder {non-binary differentiator), and interleaver, 
and a rate-1 RSC encoder (non-binary accumulator). The in- 
formation vector u is first multiplied by the transpose of P^^'', 
resulting in the intermediate vector v, with Vi ~ Ui + fl 
and = uq + /g uk^i. The vector v is then point- 
wise multiplied by the coefficient vector g'^'^' and permuted 
according to 11, obtaining a second intermediate vector v', 
which is then input to a rate-1, memory-1 time-variant RSC 
tail-biting encoder leading to a parity vector p. 



(1) ' , J-(l) 

■ 91 '"^ + fi P^-l 



Vi G [0, K - 1] 



(4) 



and with p-i = pk-i The final codeword is hence given by 
c — [u|p], and the code has length 2K symbols. The proposed 
construction brings gives a turbo code with code rate r = 
k/n ~ 1/2. We will refer to this construction as differentiate 
accumulate (DA) code. A lower rate r = 1/3 can be obtained 
by providing at the output of the encoder the intermediate 
vector V as well, i.e. by setting c = [u|p|v]. Note that, while 
for the rate 1/3 code the parity-check matrix is still given by 
^ (with proper columns permutation), for the rate 1/2 code 
([T]i represents an extended parity-check matrix, where the first 
K columns are associated with punctured symbols. The parity- 
check matrix for the rate 1/2 DA code can be obtained by 
noting that = P'^^'u^ and 



v' = vfl = u 



p(2) 



n 



The parity-check matrix is hence given by the parity-check 
matrix of a non-binary IRA code 



H 



nTp(2) p(i) 



(5) 



This parity-check matrix form is a particular case of the con- 
struction of 1271, with left and right sub-matrices characterized 
by a single cycle involving all their associated variable nodes. 



C. Cycle Graph Representation and Interleaver Design 
The codes specified by the parity-check matrices of ([T} and 



(|5]l can be conveniently described as cycle (circuit) codes 
ll27l . 1281 . as the corresponding Tanner graphs have a regular 
VN degree (i„ = 2. A graph representation for cycle codes 
can be obtained by associating a vertex with each parity-check 
equation and an edge with each codeword symbol. Considering 
the parity-check matrix ([TJ, the graph is hence given by 3K 
edges connecting 2K vertexes. 

An example is provided in Fig. |2] The graph of a rate 1/3, 
A' = 5 code (Fig. |2a)) is obtained by connecting two length- 
5 cycles according to the interleaver generated by the relative 
prime rule TT{j) = {a+p ■ j) mod K with a = 1 and p = 2. 
The graph of Fig. W[a.) turns to be the Petersen graph ll28l . ll29l . 
and hence is a (3, 5)-cage l29l . Il30l , i.e. a graph with minimal 





Fig. 2. Example of graphs for the cycle codes: (a) associated with the parity- 
check matrix of (T) for a rate 1/3, K = 5 PCCC, and (b) associated with 
the parity-check matrix of ^5) a the rate 1/2, K = 5 DA code. 



number of vertexes with degree 3 having girth 50 The graph 
associate with the parity-check matrix of (|5]l can be directly 
obtained by pruning the graph of Fig. |2a) as follows: Each 
edge connecting a vertex of the lower pentagon to a vertex of 
the upper pentagon is eliminated, and the corresponding upper 
and lower vertexes are merged together. The obtained graph 
is shown in Fig. |2b)@ 

In general, the connections between the upper and the lower 
cycles in the rate 1/3 graph define the interleaver, which may 
be selected according to rules for increasing the interleaver 
spread (turbo code perspective), or may be generated by 
filling the sub-matrix 11 of H according to girth optimization 
techniques (LDPC code perspective). The first approach has 
the inherent advantage of allowing code constructions for 
various block sizes on-the-fly by adopting efficient high-spread 
interleaver construction algorithms llSTI . ll32l FI 

III. MAP Decoding of the Component Codes 

As for the code construction, both the LDPC and the 
turbo code perspectives can be used to perform iterative 
decoding. For the former case we refer to the vast literature 
on fast Fourier transform (FFT)-based BP decoders for non- 
binary LDPC codes (see for instance ll33l ). giving decoding 
algoiithms with complexity that scales as ©(glogj q). For the 
latter, the conventional turbo decoding algorithm based on the 
BCJR algorithm ll34l applied on the trellis of the component 
codes can be simplified by FFTs as well, resulting in a 
complexity growth C'(glog2 q) as for the LDPC BP decoder 
case ITSl . We shall focus next on the symbol maximum a- 
posteriori (MAP) decoding for the component convolutional 
codes. We discuss the case of a time-variant memory-1 RSC 
code, the non-recursive convolutional code case derivation 
being similar. Each of the two RSC encoders of the proposed 

^Note that the girth of the Tanner graph associated with a cycle code is 
twice the girth of the cycle code graph, i.e. for the code based on the Petersen 
graph the Tanner graph girth is g = 2 ■ 5 = 10. 

^It is worth to note that by construction the graph of the cycle DA code is 
always given by two nested Hamiltonian cycles associated with the right and 
the left sub-matrices of the parity-check matrix (|5)- 



^ Additionally, the coefficients of g(^',f(^^ 
according to the technique introduced in (3]. 



(^',f'^' may be optimized 



PCCC scheme is in fact a time-variant memory- 1 encoder 
(i.e., the current output symbol depends only on the current 
input symbol and on the past output symbol). The time variant 
nature of the component codes is due to the multiplications 
by the coefficients gi , fi (we omitted the superscript indicating 
the branch index). The RSC encoder is fully specified by the 
relations Si = giUi + fiSi-i and pi ~ Si, where the state of 
the encoder is defined by the value stored at the input of the 
delay unit. The number of states in the code trellis corresponds 
to the field order q. 

We first consider the case where the component RSC codes 
are terminated. The a posteriori probability mass function 
(p.m.f.) vector for the symbol Ui given the channel output y is 
denoted by = [Lf (0), L^(a), . . . , if (a'"^)] with 

= Pr{ui = wjy}, G Fg. The channel observation y 
is given by y = (yo, yi, • • ■ , Yk) where each element can be 
further split as y^ = (y", y^), (j/,^) being the channel output 
corresponding to Ui (pi)u We further introduce the notation 

y[i:j] = (yi, y»+i, • • ■ , yj) (O < « < i < K). 

The computation of the a posteriori probability for the 
symbol Ui can be accomplished by evaluating 

-^"(w) = Pr{u, = uj\y] = ^ (^,_i(s)7j(s, s')|3^{s'). 

I 

The operator r"(s, s') returns the label associated with u,; for 
the trellis edge connecting the state s at time i — 1 to the state 
s' at time i, (pi-i{s) denotes the forward metric for the state 
s at time i — 1, f3i{s') is the backward metric for the state s' 
at time i, and 7i(s, s') is the transition metric between states 
s, s' at time i. We normalize the metrics such that 

ifi,{s) = Pi{S.i = s|y[o:i]}, ipo(O) = 1, 

= Pr{5, = .s|y[,+i.K]}, /3a'(0) = 1, 
7i(s,s') = Pr{S'i_i = s,5i = s'|y,} 

= PT{u,=LU,p, = iy\y'^,yP}, (6) 

with uj ~ Tl^{s,s'), V — s'. Assuming independent outputs 
UijUi^ dSll can be factored into 

7.(s, s') = Ft{u, = Lu\yn Pr{p, = i/|yf }, (7) 

7."(s,s') 7f(s') 

where 7f (s') depends on s' only since pi = Si. The for- 
ward/backward metrics can be computed recursively as 

f i{s) (x^Lpi-i{s')j,{s',s) = ^ </?,_! (s')7"(s',s)7f(s) 

s' s' 

= jf{s)J2^.,.,is'hns\s) (8) 

s' 

(3,{s) oc ^ /3,+i(s')7j+i(s,s') 

= ^A+i(,s')7."M(^>^')7f+i(5') 

= ^[/3,+i(,s')7,IVi(-^')]7r+i(5,5')- (9) 

*The vector y = (yo,yi, • ■ ■ ,yA') is composed by K + 1 elements to 
account for the additional input/output symbol required by the termination. 



Note that (O involves a convolution since s' , s are related by 
s = cjiUi + fis'. Similarly, for (|9]l s, s' are related by s' = 
5i+iUi+i + fi+is. We introduce the p.m.f. vectors 

'fit = [v»(0), fiia), • ■ • , V* (a'^"^)] 1 

A = [A(0),/3,(l), /?,(«),..., ft 

7r=[7r(0),7r(l),7r(a),.-.,7?(a^-')], 

7r = [7r(o),7r(i),7r(a),---,7r(«'^-')], 

where in the last expression (with a slight abuse of the 
notation) we re-defined 7"(w) = Pr{ui — w|y"}. In vector 
form, we can re-arrange (|8}, (|9} into 

(P^ = -ff ■ k/, ® (7^)] 

A = 7^/"+, { ■ 7f+i] ® (7r+i) } ■ (10) 

In ( [Tol l. 7ra(Q) denotes the permutation, induced by the multi- 
plication by a scalar a of a random variable with p.m.f. vector 
Q, on Q, while 7r~^(Q) denotes the inverse permutation (or 
equivalently the permutation induced by the multiplication by 
a~^). Furthermore, '•' denotes the (point-wise) multiplication 
of two vectors, and '®' denotes the convolution of the two 
vectors. The a posteriori p.m.f. vector of Ui given the channel 
output y is finally given (up to a normalization factor) by 

Lr = ^^-1 {TTf^ ® [A . 7,f]} -7," (11) 

A*" 

where p," represents the extrinsic information, fii{uj) — 
Pr{ui = Lj\y\yf}. The message update can be easily followed 
on the normal factor graph of a section of the trellis provided in 
Fig. |30 The complexity is here dominated by the convolution 
operations, and thus scales as O(g^)0 The algorithm can be 
simplified by applying the (fast) Fourier transform (FT) [|33l . 
Il36l . 1371 for finite Abelian groups on the vectors involved in 
the convolutions. Assuming extension fields with characteristic 
2, the FT reduces to the Walsh-Hadamard transform [|33l . i.e. 
given a function G F2"i, its Fourier (Walsh-Hadamard) 

transform X{i'), v G F2"i, is obtained as 

X(y)^ ^ :,(^)(_l)<-i.> 

< > being the inner product over F2 between the length- 
m binary vector representations w, of bJ,v. By employing 
FFTs, the decoding complexity is reduced to C((7log2(7). 

IV. Numerical Results 

Simulation results on the AWGN channel for codes on F256 
are presented next. In all the simulations, we adopted the BP 
decoding over the Tanner graph of the codes with a maximum 
number of iterations set to Imax ~ 200. Binary antipodal 
modulation has been considered. 

Fig. m shows the performance for a rate-compatible code 
family with input block size k = 128 bits. The mother code is 

'when decoding on tail-biting trellises , the recursion for the forward metric 
calculation shall be circularly extended (351 . 

'"Recall that q is the field order, and hence the length of the vectors involved 
in the convolutions of UP) . 
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Fig. 3. Normal factor graph for a trellis section. 



a (384, 128) code, whose parity-check matrix coefficients have 
been selected according to the method of ||3l- A lower code 
rate 1/6 has been obtained by repeating each code symbol 
twice, and by multiplying the replicas by random elements 
in F25g as for the multiplicative repeat (MR) approach of 
||2T1 . Higher code rates have been obtained in two different 
ways, i.e. (i) according to the parallel concatenation scheme, 
by periodically-puncturing parity symbols at the output of the 
two accumulators and (ii) by puncturing the VNs of type Vq 
(thus, a rate 1/2 DA code is obtained, and further higher rates 
can be achieved by puncturing symbols periodically at the 
output of the accumulator in the DA encoder). In both cases, 
symbol-wise puncturing pattern (SPP) has been applied. The 
interleaver has been designed according to a circulant version 
of the progressive edge growth (PEG) algorithm ll38l . The rate 
1/3 mother code does not show floors down to CER = 10~^, 
performing within 0.2 dB from the RGB ifTOll . Similar results 
are obtained by the lowest-rate code. For the two schemes 
with rate 1/2, the performance is still within 0.3 from the 
RGB down to CER = 10~^, with a slight advantage for the 
DA construction. The advantage is more visible for the rate 
2/3 case. Here, the PGGG performance suffers for a lack 
of steepness, which is not due to a low minimum distance 
(low-weight error patterns have not been detected), but to a 
slow decoding convergence associable with the large fraction 
of punctured symbols. For the DA case, the rate 1/2 code 
parity-check matrix of (|5]l has been used for the Tanner graph, 
and hence the higher rates have been obtained with a reduced 
fraction of punctured symbols. The same plot provides the 
performance of the (384, 128) double-binary turbo code of 
the DVB-RGS standard Q. The F256 PGGG outperforms the 
double-binary one by more that 0.7 dB at CER = 10""'. 

Fig. |5] depicts the minimum Eb/N^ required to achieve 
CER = lO"'' for several rate 1/3 parallel concatenated 
convolutional codes and rate 1/2 DA codes, with block sizes 
spanning from fc = 40 bits to fc = 1024 bits. The performances 
of rate 1/2 binary irregular protograph-based LDPG and 
accumulate repeat accumulate (ARA) codes from JSj, ||9l are 
provided too. The chart is completed by the SPB I TTI for 
the continuous-input AWGN channel. The DA codes have 
been again obtained by puncturing the Vb-type nodes of the 
PGGG graph. The interleavers have been generated on the fly 



according to |i32l . Additionally, the rate 1/3 and 1 /2 fc = 40 
bits codes associated with the cycle graphs of Fig. |2] have 
been simulated. The rate 1/3 PGGGs perform within 0.5 dB 
from the SPB all over the block sizes (with the exception of 
k = 40). For the largest (fc = 1024) block length, the gap is 
reduced to 0.3 dB. For the rate 1/2, the gap w.r.t. the SPB is 
slightly larger (0.2 dB more). The gain of the proposed non- 
binary turbo codes over the binary LDPG codes is remarkable 
1 dB or more)for the shortest block sizes. For the largest 
(fc = 1024) block length, the gain is reduced to ~ 0.3 dB. 
The performance of two short codes from ll39l are provided 
as well. The first is a (128,64,22) extended BGH code 
under maximum likelihood (ML) decoding, which achieves 
CER = lO""* at Eb/NQ ^ 3.03 dB, only - 0.3 dB away from 
the SPB with a coding gain of ^ 0.4 dB over the (128,64) 
DA code. We shall consider in the comparison that the DA 
code does not perform a complete ML decoding, and hence 
provides an error detection mechanism that may be required by 
critical application, e.g. telecommand in the up-Unk of space 
communication systems f9l|. The second code is a (600,270) 
terminated binary convolutional code with constraint length 
30. This code performs close to the (512, 256) DA code, which 
however has a slightly higher code rate (0.5 vs. 0.45) and a 
lower block size (256 vs 270 information bits). 
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Fig. 4. Performance for a rate-compatible family of turbo codes on F256, 
k = 128 bits. 



V. Conclusions 

Two novel classes of turbo codes constructed over high- 
order finite fields have been presented. The codes are derived 
from a protograph sub-ensemble of the {d^ = 2,(ic ~ 3) 
regular LDPG ensemble. One of the proposed construction is 
based on the serial concatenation of a non-binary, time-variant 
differentiator and a of non-binary, time-variant accumulator, 
and provides a highly-structured flexible encoding scheme for 
{dy = 2,4 = 4) LDPG ensembles. Symbol MAP decoding 
of the component codes has been illustrated, together with its 
FFT-based simplification. The proposed codes allow efficient 
decoding either as LDPG or as turbo codes. Remarkable 
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Fig. 5. Required E^/Nq for achieving CER = 10 * for various codes witli 
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gains 1 dB) w.r.t. binary LDPC/turbo codes have been 
demonstrated in the moderate-short block regimes. 
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